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The numerically exact diagonalization study on small clusters of the t-J model with second- 
and third-neighbor hopping terms shows that a novel spin liquid state is realized around a doped 
hole with momentum k=(7r,0) and energy ~2J compared with that with (n/2,n/2) in insulating 
cuprates, where the spin and charge degrees of freedom are approximately decoupled. Our finding 
implies that the excitations in the insulating cuprates are mapped onto the d-wave resonating 
valence bond state. 
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The pairing symmetry of high-T c superconductors has 
been established to be of d x i_ %j -i type, whose gap is 
maximum in the direction of k=(7T,0). Such a d-wave 
gap has been observed not only in the superconducting 
state but also in the normal state for the underdoped 
superconductors by angle-resolved chotoemission spec- 
troscopy (ARPES) experimentsBErtf The normal-state 
gap is called a pseudogap. In the overdoped cuprates, 
the pseudogap almost disappears and a flat band with 
sharp spectral weight emerges around the (7r,0) point. 
Understanding of the electronic states at (n,0) is, there- 
fore, critical in the field of cuprate superconductivity.. n 

In the insulating cuprates, ARPES experimental Er 
have revealed the characteristic excitations below the 
charge transfer gap with minimum binding energy at 
k=(7r/2,7r/2): While the spectrum near (ir/2,w/2) con- 
sists of a sharp peak, the (7T,0) spectrum is very broad. 
From the numerically exact diagonalization study for a 
t-J model with second- and third-neighbor hoopings t' 
and t", respectively, i.e., the t-t'-t"-J model,© it has 
been found that the quasi-particle (QP) weight at (n,0) 
is remarkably reduced through weakening of the anti- 
ferromagnetic (AF) spin correlation induced by t' and 
t". The resulting broad (n,0) spectrum in the insulator 
continuously evolves into the bijoad one at underdoping 
to the sharp one at overdoping.B Although the numer- 
ical results have demonstrated a reasonable agreement 
with ARPES data, the physical origin of the broadness 
of the (7r,0) spectrum remains to be elucidated. Very 
recently, Ronning et alu pursued this issue, perform- 
ing ARPES experiments on Ca2CuC>2Cl2, a parent com- 
pound of high-T c superconductors. The experimental 
data show a d-wave-like dispersion along the (tt/2,tt/2)- 
(•7r,0) line, implying that the pseudogap is remnant of the 
G?-wave gap in the insulator. 
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In this Letter, we examine the dynamics of a hole in 
the insulator and the spin configuration around it. With 
realistic values of t' and t", we find a dispersion consis- 
tent with that-predicted by the resonating valence bond 
(RVB) theory© with d-wave pairingJj'BEiJ where the in- 
coherent motion of the hole plays a crucial role. These 
t' and t" also lead to a novel spin liquid state around a 
hole with k= (ir,0) and excitation energy of ~2 J, where 
the hole causes an antiphase in the AF spin background. 
Correspondingly, the spin and charge degrees of freedom 
are approximately decoupled there, as seen from the spin 
and charge correlation functions. In contrast to the (7r,0) 
state, the low-energy states near (ir/2,ir/2) exhibit AF 
behavior. The implication of the spin liquid state is dis- 
cussed in connection with the d-wave RVB state as well 
as the pseudogap. 

We first clarify the importance of the incoherent mo- 
tion of a hole by comparing the calculations of the single- 
hole dispersion based_ao. a self-consistent Born ap- 
proximation (SCBA)BBB) for the t-t'-t"-J model and 
(ii) a spin-density wave (SDW) mean-field approximation 
for a Hubbard model containing the long-range hopping 
terms (t-t'-t" Hubbard). In SCBA, the self-energy of a 
hole is determined by taking into account the coupling 
to the spin background. This scheme explains consis- 
tently the reduction of the dispersion from t to J and 
the incoherent part of the spectral function extending 
over the energy ~ t. The SCBA results are in good 
agreement with those-iif the numerically exact diagonal- 
ization (ED) methodS In contrast to SCBA, the SDW 
approximation gives only a coherent motion of a hole. 

The t-t'-t" -J model is given by 

H = J E Si'Sj-t E 4r C jV 

E 4*Cj<r-t" E C .W+H.C.,(1) 
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Fig. 2. (a) The dependence of the energy difference i?(7r,0) — 
E(tt/2,tt/2) (in units of /=0.35 eV) on the values of /' and t" . 
The variable a represents a scaling factor of t' and t" as t' (t") = 
crf (t£,'), being ^=-0.12 eV and t"=0.08 eV. Solid line: SDW 
result of the t-t'-t" Hubbard model given by 4a (2t'^ - /[>) . Solid 
circles: SCBA results of the t-t'-t"- J model on a 16x16 lattice 
with J/t=0A. (b) (E(k) -B(7r/2,7r/2)) versus J/t in the case 
of a=l. The solid line represents 2.2J/t known as the band 
width from (7r/2,7r/2) to (0,0). 



Fig. 1. (a) Contour plot of SCBA single- hole dispersion of the ( 
t'-t"-J model on a 16x16 lattice with t=0.35 eV, t'=-0.12 c\ 
t"=0.08 eV and J/t=0.4. The numbers in the plot are energ 
values in units of eV relative to the (n-/2,7r/2) point, (b) Th 
same as (a) but for SDW mean-field approximation results of 
the half-filled t-t'-t" Hubbard model with the same hopping pa- 
rameters and U /t=10. (c) d-wave plot of the dispersions along 
(7r/2,7r/2) to (tt,0). Dashed line: SDW result of the t-t'-t" Hub- 
bard model. Solid circles: SCBA results of the t-t'-t"-J model. 
Open squares: ARPES data of Ca2Cu02Cl2 taken from ref. 6. 



where the summations («,i) lst , j) 2n d anc ^ (*>j)3rd run 
over first-, second- and third-nearest-neighbor pairs, re- 
spectively. No double occupancy is allowed, and the 
rest of, the notation is standard. From the previous 
study IP it has been known that £=0.35 eV, t'=— 0.12 eV, 
and i"=0.08 eV for Sr 2 Cu0 2 Cl 2 and Bi 2 Sr 2 CaCu 2 8+< 5. 
These values reveal consistent results of the Fermi surface 
topology with experimental data. The SDW dispersion 
inJic AF state of the t-t'-t" Hubbard model is given 
bye 

E±=e k ±E[-», (2) 

with 

e^. = —At' cos k x cos k y — 2t" (cos 2k x + cos 2k y ) , (3) 

= \J& 2 (cos k x + cos k y f + (Umf. (4) 

The chemical potential \x and the magnetization m are 
determined by solving a set of self-consistent equations. 
We use the same values of t, t' and t" as those for the 
t-t'-t" -J model, and U/t=10 (U is the on-site Coulomb 
interaction) . 

Figures 1(a) and 1(b) show the contour plots of the 
single-hole QP dispersion obtained by the SCBA for the 
t-t'-t" -J model on a 16x16 lattice and by the SDW 
approximation for the t-t'-t" Hubbard model, respec- 
tively. The SCBA contour is isotropic and thus consis- 



tent with the ARPES dispersion in Ca 2 Cu0 2 Cl 2 that 
is isotropic around (ir/2,ir/2), while the SDW one is 
not. In addition, the width of the SDW dispersion 
is remarkably wider than that for the SCBA results. 
This is clearly seen when we plot the dispersion along 
(tt/2,tt/2) to (7r,0). In Fig. 1(c), the energies relative 
to the (7r/2,7r/2) point, E(k) - E(tt/2, 7t/2), are plotted 
against \cosk x — cosk y \ /2, together with the data from 
Ca 2 Cu0 2 Cl 2 which are well fitted by a straight lincB 
What is the origin and implication of the difference be- 
tween the two results? Figure 2(a) shows the dependence 
of the energy difference E(ir, 0) — E(ir/2, tt/2) on the val- 
ues of t' and t" . The variable a represents a scaling factor 
of t' and t" as t' = at' and t" = at%, being ^=-0.12 eV 
and £q=0.08 eV. The SDW result is given by a function 
4a (2t# - t' ) [see Eq. (2)], independent of U. The result 
of the t-t'-t" -J model deviates remarkably from the SDW 
line and has a tendency to saturate around a ~ 1. In 
addition to the slow increase of the energy with increas- 
ing a, the wave-function renormalization Z (ref. 16) at 
(7r,0) is dramatically suppressed from Z=0.87 (a=0) to 
0.24 (a=l). This is again in contrast to the SDW result, 
where Z is always one. At a ~ 1, the energy difference is 
approximately 0.8£=2 J. This number is of crucial impor- 
tance because the incoherent spectrum of a hole starts 
from the energy position of about 2 J above the spectrum 
of QP at (7r/2,7r/2) (see insets in Fig. 4). Although the 
energy of the QP spectrum at (tt,0) is almost the same 
as that at (tt/2,tt/2) in the t-J model, the former in- 
creases with increasing a and the spectrum merges into 
the incoherent spectra for a~l. As a result, the differ- 
ence between E(n,0) and E(tt/2,tt/2) saturates and Z 
becomes small. 

The above results indicate that when realistic values 
of t' and t" (a— I) are employed in the t-t'-t" -J model, 
the dispersion from (ir/2,ir/2) to (tv,0) is sensitive to the 
value of J. We now plot, in Fig. 2(b), the energy differ- 
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Fig. 3. (a) Spin correlation Cggi (k) of the t-t'-t"-J model on a 2C 
site cluster with one hole. t=0.35 eV and J/t=0A. The variabl 
a represents a scaling factor of t' and t" as t'(t") = at' (t^ ,, 
being t' =-Q.l2 eV and t"=0.08 eV. k=(7r,0). (b) The same as 
(a) but for k=(0,0). (c) The same as (a) but for k=(27r/5,7r/5). 
The labels in (a), (b) and (c) denote configurations around a 
hole shown in the inset of (a), (d) The staggered spin correlation 
D^(r) at several momenta for a=0 and 1 of the 20-site t-t'-t"-J 
cluster as a function of |r|. The result for the Heisenberg model 
is also shown for comparison. 



ence E(ir, 0) — E(ir/2, ir/2) as a function of J/t, together 
with the difference between E(0, 0) and E(n/2, tt/2) that 
has been known to be scaled as ~ 2. 2 J. if We find that 
E(ir, 0) — E(tt/2,tt/2) has the same J dependence as 
£(0,0) — E(n/2,n/2), at least, in the realistic range 
J/i=0.2~0.6. Therefore, the dispersion of the hole is 
governed by the spin degree of freedom. This explains 
why the dispersion is isotropic around (ir/2,ir/2). Such 
an isotropic dispersion whose width is controlled by J 
is similar n t£L that from the flux phase or d-w&ve RVB 
pictureBBe 

We now examine the spin and charge characteristics 
of the QP state with k=(7r,0). Figure 3(a) shows the 
spin correlation around a hole defined as Cgg' (k) = 
£V (</>Qp(k) |n^S i+ a • S i+5 < | c/> QP (k)), as a function of a. 
Here, 5 and 5' denote two sites around the hole follow- 
ing the labeling convention shown in the inset, n\ is the 
hole-number operator at site i, and </>Qp(k) represents 
the wave function of the QP state with momentum k. 
The numerically exact diagonalization method is used 
for a 20-site cluster. For a=0, i.e., the t-J model, the 
spin correlations in the a and b configurations [see the 
inset of Fig. 3(a)] are positive, i.e., ferromagnetic (FM), 
while those for d and e are AFJB> With increasing a, 
the correlation of a changes from FM to AF and satu- 
rates near a ~ 1. The b configuration shows a similar 
but slow change. The configuration d exhibits the oppo- 
site behavior, reducing AF correlation with increasing a, 
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Fig. 4. A(k, uj) at half-filling, and 5(q, a;) and 7V(q, ui) with q=k- 
ko at one-hole doping for a 20-site t-t'-t"-J model. J/t=0A. 
k=(7r,0) and ko = (37r/5,— 7r/5). The energy w in j4(k, u>) is mea- 
sured from the QP peak at ko- (a) t-J and (b) t-t'-t"-J with real- 
istic values oft' and t" (a=l). The height of the bars represents 
the spectral intensity. The insets in (a) and (b) show A(k, w) 
with full-energy scale. The delta functions (vertical bars) are 
broadened by a Lorentzian with a width of 0.3t (solid lines). 



which again nearly saturates around a=l. These behav- 
iors from a=0 to ~1 imply that the spin background of 
the (7r,0) state changes from a Neel-like state to a spin 
liquid state in which spins have no tendency to order. It 
is interesting that the spin liquid state is accompanied by 
an antiphase of spins around the hole, as evidenced by 
AF correlation of the a configuration for Since the 

correlation of b is very small there, the antiphase of spins 
is one-dimensional (ID). This situation is similar to the 
case of a doped ID. insulator in which the spin-charge 
separation occurs.li3> In contrast to the (nfl) state, QP 
states with different momenta do not show such remark- 
able changes as a function of a, as shown in Figs. 3(b) 
and 3(c). We note that the AF correlation of the a config- 
uration for k=(27r/5,7r/5) is not related to the spin liquid 
state [see also Fig. 3(d)] but due to a twisted characteiurf 
the spin correlation seen in the one-hole ground state. t3 

The staggered spin correlation D u(r) = J2i Pi( r ) 
(</)Qp(k) \ S i+r ■ Si | </>Qp(k)) /N s also supports the con- 
cept of the spin liquid state with (tt,0). Here, r is the 
vector connecting two sites, N s is the total number of 
sites, and i^(r)=±l depending on whether the two sites 
are on the same sublattice. Figure 3(d) shows D-^(r) 
at several momenta for a=0 and 1, together with the 
spin correlation of the Heisenberg model which shows 
AF long-range order in the thermodynamic limit. Ex- 
cept for the case of a—1 at (tt,0) (filled circles), D-^(r) 
shows |r| dependence similar to that of the Heisenberg 
model, indicating the presence of long-range order. The 
(7r,0) state for a=l, however, shows a rapid decrease of 
the correlation with increasing |r|. This is again consis- 
tent with the spin liquid concept. 

The dynamical properties of spin and charge degrees of 
freedom also provide us with useful information on the 
QP state with (ir,0). We calculate the dynamical spin 
and charge correlation functions, S^q, u>) and iV(q, w),E3' 
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on the 20-site cluster with one hole, the Hilbert space of 
which is equivalent to that for the final state of the single- 
hole spectral function A(k, oS). The momentum transfer 
is given by q=k— ko, where ko is the ground-state mo- 
mentum of the one-hole system. In Fig. 4, A(k, oj) with 
k=(7T,0) is compared with_5(q, w) and N(q,w), where 
we take k =(37r/5,— 7r/5).E$ In the t-J model, both spin 
and charg£_components are involved in the QP state at 
cj/t=0.08,EiP while only spin component remains in the 
lowest-energy state at u/t—QAl in the t-t'-t"-J model 
as seen in Fig. 4(b). This is because of the separation 
of the spin and charge degrees of freedom. However, 
the separation is incomplete unlike the ID case of the 
t-J modelM We note that the results of the t-t'-t"-J 
model are smoothly connected to those of the t-J model, 
showing no separation, as seen in Fig. 4(a). Much of 
the weight of the charge excitation shifts to a higher en- 
ergy region where the spectral function has large weight 
(u/t ~7) [see the inset of Fig. 4(b)]. For k=(27r/5,7r/5), 
which is close to (tt/2,it/2), the QP involves both spin 
and cJaarge components independent of a (not shown 
here)B) 

We have shown that the spin background in the (tt,0) 
state in the insulator behaves like a spin liquid state. At 
the same time, the spin and charge degrees of freedom 
are nearly separated in the stateO The almost linear 
behavior of the dispersion near (ir,0) in the d-wave plot 
shown in Fig. 1(c) indicates that the spin liquid state has 
a (i-wave gap. All of these facts tempt us to identify the 
state with the d-wave RVB state discussed in refs. 8—10. 
However, we have to distinguish the present spin liq- 
uid state with the RVB state: the former is seen in the 
one-hole state with k=(7r,0) and the excitation energy of 
~2J, while the d-wave RVB theory predicts a spin liq- 
uid state as well as the spin-charge separation which are 
independent of the momentum of a doped hole. In fact, 
the present results reveal that hole states near (n/2,ir/2), 
i.e., the low-energy states, are not effective to destroy 
the AF long-range order [see Fig. 3(d)]. However, the 
states near (ir,0), i.e., the high-energy (~ 2J) states, are 
in good agreement with the excited states predicted by 
the d-wave RVB theory. Therefore, the spin background 
changes from the AF state in the low-energy— region to 
the spin liquid one in the high-energy region^ The d- 
wave RVB theory might thus be regarded as an effective 
theory to describe the states in the high-energy region. 

With hole doping, the energy of the (7r,0) r state shifts 
to the Fermi level, maintaining the broadnessB that is an 
indication of the spin liquid state. Since the pseudogap 
with the order of J (high-energy pseudogap) is deter- 
mined by the position of the the broad spectrum,© the 
spin liquid concept is the most probable explanation of 
the origin, of the high-energy pseudogap in underdoped 
cuprates.tf 

In summary, we have examined the dynamics of a 
doped hole and spin correlation around it in insulating 
cuprates by using the t-t'-t"-J model. In contrast to the 
AF state near (tt/2,it/2), we have found a novel spin liq- 
uid state around the hole with k=(7r,0) and excitation 
energy of ~ 2 J, where the spin and charge degrees of 
freedom are approximately decoupled. Our findings im- 



ply that the excitations in the insulating cuprates are 
mapped onto the d-wave RVB state. 
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